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, PACS. 83.80.Va - Elastomeric polymers in Rheology. 

PACS. 46.70.Lk - Structures other than beams, plates etc. in Continuum mechanics of solids. 
PACS. 81.05.Lg - Polymers and plastics in Material Science. 
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Abstract. - After thrown onto a solid substrate, very soft spherical gels bounce repeatedly. 
Separate rheological measurements suggest that these balls can be treated as nearly elastic. 
The Hertz contact deformation expected in the static (elastic) limit was observed only at very 
small impact velocities. For larger velocities, the gel ball deformed into flattened forms like a 
pancake. We measured the size of the gel balls at the maximal deformation and the contact 
time as a function of velocities for the samples different in the original spherical radius and the 

S Young modulus. The experimental results revealed a number of scaling relations. To interpret 
these relations, we developed scaling arguments to propose a physical picture. 
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Introduction. - Impact of spherical objects is concerned with a wide range of subjects 
from daily life to solid mechanics [TJEJG]. Accordingly, it has been studied for a long time in 
different contexts and still continues to be an important research subject. The Hertz contact 
theory 0] is a classic example in the static limit within the linear elastic theory. For larger 
impact velocities, internal vibrations [S] and dissipative processes such as solid viscosity 0, 
plastic deformation should be considered. Impact of solid spheres in a granular medium 
has been found to be important fundamental processes in granular flow, and has attracted 
a considerable attention jTHSUHl ■ Impacts of balls of soft matter could also be an interesting- 
problem: unique phenomena are expected due to a variety of constitutive equations and due 
to large deformations. For example, impact or dynamics of liquid balls has attracted a wide 
audience [I1II33E1II3E3E1E1E7I . 

In this paper, we study another example from soft matter: impacts of soft gel balls verti- 
cally impinging onto a solid substrate. There the Hertz contact deformation is observed only 
for small velocities. For larger velocities, gel balls are flattened globally during the impact on 
the substrate, i.e., expanded into the lateral direction. We present experimental data on the 
lateral dimension of balls at the maximal deformation and the contact time as a function of 
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sample 


Water 


AA 


BIS 


E(Pa) 


BIS4 


lOOcc 


lOg 


0.04g 


1.24x10* 


BISIO 


lOOcc 


lOg 


o.ig 


2.71 xlO 4 


BIS15 


lOOcc 


lOg 


0.15g 


3.87xl0 4 


BIS20 


lOOcc 


lOg 


0.20g 


4.56xl0 4 



Table I - The composition of four samples of acrylamide gel and their Young moduli E. 



impact velocities to show the existence of a number of scaling relations. These relations are 
interpreted via a naive physical picture. 

Sample gels. - We used a series of acrylamide gels with the same polymer concentration 
but with different cross-link densities; acrylamide monomer (AA, M w =71.08) constitutes 
sub-chains while methylenebisacrylamide (BIS, M w =154.17) cross-links. The amount of each 
reagent for preparing acrylamide gels is shown in Table 1. Ammonium persulphate (lwt % 
of AA) and tetramethylethlylcncdiamine (0.25vol % of water) were added to initiate and to 
accelerate the radical polymerization of A A and BIS. The pre-gel solutions were sealed into 
a spherical mold (R =14mm or 31.5mm) consisting of two hemispherical shells. Gelation 
reaction continued for 24 hours at 30 °C. 

Rheological characterization was carried out on cylindrical gels with a rheometer (RE- 
OGEL, UBM Co.) in an oscillatory compression mode. As a result, we found that the 
mechanical responses of our samples can be all regarded as nearly elastic. The values of tan 6 
are less than 0.02 and the real part of the complex modulus E 1 is almost constant in the range 
of strain frequency between 0.1Hz and 100Hz for all samples except BIS4 (Here, 5 is the phase 
difference between the stress and the strain). Even for this softest sample, tan 5 increases 
by a fairly small amount with / (for example, tan 5 ~ 0.02 at / = 0.01Hz, tan 6 ~ 0.05 at 
/ = 50Hz and tan 5 ~ 0.1 at / = 100Hz). Tabled also shows the modulus E' at / = 10Hz 
that is regarded as the static Young's modulus E in the following. 

Impact experiment. - Gel balls freely fall on a fixed aluminum plate of 20mm thickness. 
Before the impact, the gel ball is pinned at a height L from the plate by a tube sucking air 
weakly. The gel ball begins to fall by switching-off the sucking, and then impinging onto the 
aluminium substrate. We coated the surface of balls in white with aluminum oxide powder 
to avoid sticking of the gel balls to the substrate (this is especially important for low impact 
velocities). The impact velocity V can be determined from the relation V = y/2gL. The 
impact processes are recorded by a high-speed CCD video camera (Motion Coder Analyzer 
SR: Kodak Co.) with recording rates of 1000 FPS (samples of R = 31.5mm) or 2000 FPS 
(samples of R = 14mm). 

Experimental Results. - We found experimentally that modes of deformation are rather 
different depending on the velocities. The form of maximal deformations can be categorized in 
three classes: (1) the Hertz type for very small impact velocities (Fig. QJi), (2) quasi-ellipsoid 
for intermediate velocities (Fig. HJ)), and pancake for large velocities (Fig. QJ;). 

For small impact velocities the Hertz-type deformation is theoretically expected: in the 
static limit, because of the nearly pure elasticity of gel balls, the localized Hertz deformation 
should be observed. However, this theory should break down for large impact velocities where 
a ball deforms non-locally. Indeed, the Hertz regime can be observed only for very small 
velocities: the impact in Fig. was achieved by dropping the gel ball from a low height 
L less than 1cm (without using the sucking system) where a precise determination of V is 
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Fig. 1 - Maximal deformations of gel balls (R — 31.5mm). a) Hertz regime for very low impact 
velocities (BIS10, V < 0.4m/s). b) quasi-ellipsoid regime for intermediate velocities (BIS4, V = 
2.0m/s). c) pancake regime for high velocities (BIS4, V = 7.0m/s). 



difficult (V < 0.4m/s) (We note here that all the data points used in the following plots are 
limited to the case where velocity determination can be done rather precisely). For moderate 
impact velocities deformation is no longer localized around the contact area and the ball takes 
a flattened shape which is quasi-ellipsoidal (Fig. ^>). For large impact velocities the ball is 
strongly deformed to take a pancake shape (Fig. H 3 ). 

We measured the horizontal radius R + X (see Fig. ^| at the maximal deformation and the 
contact time at various conditions. As shown below, we found that these experimental data 
can be well characterized by the length scale R, the sound velocity V c and the corresponding 
time scale r c : 

V c - (1) 



t c = R/V , (2) 

where p is the density of gel balls (p ~ 1.05 x 10 3 kg/m 3 ). Characteristic scales V c and r c come 
out naturally from theoretical considerations presented below. 
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Fig. 2 - Maximal deformation X/R as a function of impact velocities V/V c for (a) smaller samples 
(R =14mm) and (b) larger samples (R —31. 5mm). Two dashed lines in the plots corresponds to lines 
with slopes 1 and 8/5, respectively. 

Fig. |21 shows the maximal deformation X/R as a function of reduced velocities V/V c . As 
mentioned above, owing to the characteristic scales R and V c , the data from the four samples 
for each size (R = 14mm or 31.5mm) collapse well onto a single behavior which can be divided 
in two regimes as indicated by two dashed lines with different slopes (8/5 and 1). For large 
velocities (V > V c ) it scales as X ~ V while for small velocities (V > V c ) X ~ V n with n 
being a certain value larger than the unity. A theoretical value of the exponent, n = 8/5, for 
the small velocity region seems consistent with the data (see below). 
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(b). 



3 - Contact r as a function of velocities V for gel balls of radius R —14mm (a) and R =31. 5mm 



Fig. [3]shows a contact time r as a function of velocities. The two plots collapse well onto a 
single master curve in Fig. thanks to the characteristic scales r c and V c as mentioned above. 
The contact time r decreases with increase in velocities, following a scaling law (r ~ V -1 / 5 ) 
and seems to approach a constant value (Figs. |31and0J). The plateau value is of the order 
of millisecond. This value increases with the size R and decreases with the modulus E. 
Accordingly, the data from different sizes and moduli collapse by the characteristic scale r c . 

As seen above, both X and r can be divided in two regimes by a characteristic velocity of 
the order of m/s. This velocity scale increases with a modulus E (this can be confirmed due 
to the data collapse via renormalization of velocity by V c ). 

Theoretical considerations. - As mentioned above and as seen from the snapshots at the 
maximal deformation, modes of deformation seem to change with increase in impact velocity: 
from (1) Hertz type, to (2) ellipsoid, and then to (3) pancake (see Fig^l. We shall develop a 
theory based on these observed forms of deformation. We note in advance that experimental 
results can be understood by combining only the two limiting regimes (Hertz and pancake). 
However, we still present the result based on the ellipsoid approximation as reference, which 
is consistent with the experimental results in the turn-over regime. We start from scaling 
arguments for the Hertz regime only to reproduce the essential properties of the well-known 
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Fig. 4 - Reduced plot of a contact time r versus impact velocity V. The right plot suggests a scaling 
relation r ~ V~ 1/5 . 
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complete analytical solution. The arguments are then extended to the other regimes. 

In the Hertz regime at small impact velocities, the elastic energy stored at the maximal 
deformation can be estimated dimensionally as Ujj ~ E [Y/a) 2 a 3 with E being the elastic 
modulus of the ball (Y is defined as in Fig^and a is the radius of the contact area); a strain 
field of the order of Y/a is localized in a volume of the order of a 3 (the deformation field 
relaxes out at a distance of the order of a from the contact surface because the strain field is 
governed by the Laplace equation). Due to the geometrical relation Y ~ a 2 / R with R being 
the radius of the ball, the energy reduces to the well-known form: 

rb/2 



U H ~ EVRY J/2 (3) 

Assumption of the energy conservation between the initial time and the maximally-deformed 
moment M V 2 ~ Uh (M ~ pR 3 is the mass of the ball) leads us to an estimate for the contact 
time r ~ Y/V, which results in the well-known relation, 

t/tc-OTO-VS (4) 

where the characteristic scales V c and t c given in Eqs. (JIJ and have naturally come out 
as announced. This implies a moderate increase of the contact time with decrease in velocity. 
To obtain a non-trivial X — V relation, we assume that in our " Hertz regime" the shape is a 
sphere cut out by a plane (see below) and require the incompressibility condition: the part of 
the sphere cut out by the substrate is compensated by an increase in radius, i.e., 

a 2 Y ~ R 2 X. (5) 

Combined with the previous relation of the energy conservation, we have the relation between 
impact velocity and deformation in the horizontal direction as mentioned above: 

X/R ~ {V/V c f 5 (6) 

In the ellipsoid regime at intermediate velocities, the maximal elastic energy is given by 
Ue ~ E (Y/R) 2 R 3 ; a strain field of the order of Y/R ~ X/R is distributed within the whole 
volume of the order of i? 3 . Note here the relation, Y = 2X, which expresses the condition of 
volume conservation for small deformation (X/R < 1). Thus, we have 

U E ~ ERY 2 (7) 

The energy conservation, MV 2 ~ Ue, allows us to obtain an estimate for the contact time 
and a velocity-radius relation, 

T ~ T C (8) 

X/R ~ V/V c (9) 

In the pancake regime at high velocities, where X > R > Y, the maximal energy becomes 
Up ~ E ((R + X)/R) 2 R 3 (the ideal rubber deformation energy), or a linear spring energy 

Up - ERX 2 

Noting that, in this case, the contact time might well be estimated not by Y/V but by X/V 
due to an analogy with a spring system (physically, we can imagine that after a strong impact 
the vertical velocity is immediately redirected toward the horizontal direction) , we obtain the 
same scaling relations with the ellipsoid case (Eqs. © and ©)■ 
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The transition between the Hertz and the ellipsoid regimes is given by the condition, 
Ue ~ Uh, which implies Y < R. In terms of velocity this is expressed as V ~ V c (compare 
Eqs. J3J and JSJ). On the other hand, the transition between the ellipsoid and the pancake 
regimes should be marked by X > Y, or V > V c . Since the transition from Hertz to ellipsoid 
regimes and that from ellipsoid to pancake regimes are predicted to occur at the same velocity 
at the level of scaling laws, the ellipsoid regime is not expected to manifest itself clearly. 

We summarize the above theoretical predictions which agree with experimental observa- 
tions: (1) X scales as X/R ~ (V/V c ) s ^ 5 for small velocities but as X/R ~ V/V c for large 
velocities. (2) The contact time r is constant for large impact velocities but below a certain 
velocity it deviates from the plateau value and increases with velocity decrease following a 
scaling law, r ~ V^ 1 ^ 5 . Typical plateau value r c for our samples is of the order of a few ms. It 
increases with R while decreases with E. (3) Both X and r can be divided in the two regimes 
by a transition velocity V c . Typical value V c is of the order of a few m/s and it increases with 
E. (4) Data should be well characterized by reducing size, time, and velocity variables by 
scales R, r c , and V c , respectively. 

In our Hertz regime, we assumed that the shape is a sphere cut out by a plane and, in 
addition, we required Eq. JSJ, which is natural at least in the case of water drops where the 
surface energy is dominant |13U18j . As a result, we have a global strain X/R distributed over 
the whole volume R 3 , in addition to the original local distribution Y/a over a 3 . The former 
extra energy is far smaller than the latter with a small ratio ~ (a/i?) 3 ; this correction does 
not change the results of our previous analysis. 

Conclusion. - Experimentally, the scaling relation X ~ V for large velocities is shown 
and this relation seems to change into another scaling relation with a stronger power below 
V ~ V c - Another scaling relation r ~ V^ 1 ^ 5 for small velocities is also shown while, for larger 
velocities (V > V c ), t seems to approach a plateau value. These behaviors of r and X can 
be explained by the theory which starts from experimentally observed shapes of balls at the 
maximal deformation and employes the energy conservation. This theory also leads to the 
prediction X ~ I/ 8 / 5 for small velocities, which seems consistent with experimental data. 

Discussion. - Experimentally it is difficult to observe a wide plateau region of the con- 
tact time predicted by the pancake form and only an asymptotic behavior (or the ellipsoid 
behavior) to this limiting regime is observed. This is because of (1) it is practically difficult to 
achieve such high impact velocities by the present experimental setup and (2) it is inherently 
impossible to exceed a certain high impact velocity above which the impact causes irreversible 
damages to gel balls. 

The effect of gravity can be another source of increase in contact time with decrease in 
impact velocity as might be the case in some bouncing water drops In the present case, 
however, this possibility seems to be excluded; the ratio MgY/Un ~ pgR(R/Y) 3 ^ 2 / E suggests 
that the gravitational energy MgY becomes important only for Y < 0.04_R, which is outside 
of our experimental region. Furthermore, if we define the length l g by pgl g ~ E, which is 
a counterpart of the capillary length, this is about lm and is well beyond the characteristic 
length scale R (of the order of cm). This also suggests a weak gravity effect. These arguments, 
in turn, suggest that the gravity may play a role for very small deformations, which are not 
studied here. 

The effect of viscosity is possibly more important and will be discussed in detail elsewhere. 
In fact, if we closely look at Fig. 0]we find that the behavior of the softest sample with the 
highest dissipation seems to deviate slightly from the others (if we removed these points from 
these plots, the data collapse would be much better). This suggests a possibility that viscous 
effects become important already for the softest sample. 
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The characteristic velocity V c that emerged as the result of the static energy evaluation 
turned out to coincide with the velocity of the shear wave propagation, at the level of scaling 
laws. To understand this meaning, a systematic study on the static large deformations of soft 
gel balls would be helpful, and thus such measurements are now under study. 

Acknowledgements. - Y. T. thanks Yuji Kozono and Yoshi Miyamoto for their instruction 
in the rheological measurements. Y. T. also thanks Takanobu Sato for his cooperation at the 
early stage of this study. Y. T. and Y. Y. thank Mitsugu Matsushita for his encouragement. 
K. O. benefited a lot from David Quere through useful discussions. K. O. is very grateful to 
P.-G. de Gennes for his third stay in Paris financially supported by College de France. K. 
O. also appreciates Christophe Clanet and Frederic Chevy for discussions. This work is also 
supported by an internal grant of Ochanomizu University. 

REFERENCES 

[1] JOHNSON K. J., Contact Mechanics (Cambridge University Press, Cambridge, 1960). 

[2] GOLDSMITH W. J., Impact (Edward Arnold Pub., London, 1960). 

[3] STRONGE W. J., Impact mechanics (Cambridge University Press, London, 2000). 

[4] HERTZ H., J. Reine Angew. Math. 92 (1882) 156. 

[5] GERL F. and ZIPPELIUS A., Phys. Rev. E, 59 (1999) 2361. 

[6] KUWABARA G. and KONO K., Jpn. J. Appl. Phys., 26 (1987) 1230. 

[7] LABOUS L., ROSATO A. D. and DAVE R. N., Phys. Rev. E, 56 (1997) 5717. 

[8] MORGADO W. A. M. and OPPENHEIM I., Phys. Rev. E, 55 (1997) 1940. 

[9] HAYAKAWA H. and KUNINAKA H., Chem. Eng. ScL, 157 (2002) 239 and references therein. 
[10] BERGERON V. and QUERE D., Phys. World, 14 (2001) 27. 
[11] RICHARD D., CLANET C. and QUERE D., Nature (London), 417 (2002) 811. 
[12] RICHARD D. and QUERE D., Europhys.Lett, 50, 769(2000). 
[13] MAHADEVAN L. and POMEAU Y. Phys. Fluids, 11 (1999) 2449. 

[14] BERGERON V., BONN D., MARTIN J. Y. and VOVELLE L., Nature (London) 405 (2000) 
772. 

[15] PEREZ M., BRECHET Y., SALVO L., PAPOULAR M., and SUERY M., Europhys. Lett. 47 
(1999) 189. 

[16] AUSSILLOUS P. and QUERE D., Nature (London), 411 (2001) 924. 
[17] MAHADEVAN L., Nature (London), 411 (2001) 895. 

[18] De Gennes P.-G, Brochard-Wyart F., and QUERE D., Gouttes, bulles, perles et ondes (Belin, 
Paris, 2002). 

[19] OKUMURA K., CLANET C, RICHARD D. and QUERE D., Europhys. Lett. 62 (2003) 237. 
In the case of water drop, we cannot consider the Hertz type deformation as the source of the 
contact time increase with decrease in velocities because such deformation mode results in the 
constant contact time (note that there are no elastic contributions and the surface increase 
scales as X 2 ). 



